Kondo screening cloud in the single-impurity Anderson model: A DMRG study 
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A magnetic moment in a metal or in a quantum dot is, at low temperatures, screened by the 
conduction electrons through the mechanism of the Kondo effect. This gives rise to spin-spin corre- 
lations between the magnetic moment and the conduction electrons, which can have a substantial 
spatial extension. We study this phenomenon, the so-called Kondo cloud, by means of the density 
matrix renormalization group method for the case of the single-impurity Anderson model. We focus 
on the question whether the Kondo screening length, typically assumed to be proportional to the 
inverse Kondo temperature, can be extracted from the spin-spin correlations. For several mecha- 
nisms - the gate potential and a magnetic field - which destroy the Kondo effect, we investigate the 
behavior of the screening cloud induced by these perturbations. 

PACS numbers: 78.20.Bh, 02.70.-c, 72.15.Qm, 75.20.Hr 



I. INTRODUCTION 

The Kondo efFect^i a well-known feature of magnetic 
impurity systems, has seen a tremendous renewed inter- 
est due to the realization of quantum dots and nanoscale 
systems 1^ The existence of Kondo correlations at low 
temperatures T has been firmly established in numer- 
ous experiments on quantum dots,— molecules,— and car- 
bon nanotubes.'^ The interaction of an impurity spin 
with itinerant electrons, causing the Kondo effect, man- 
ifests itself in spatially extended spin-spin correlations 
- the Kondo screening cloud. These correlations have 
been extensively studied in theory^i^ i^i^i-^^d^i^^d^'-^^'-^^ 
and many proposals for experimentally measuring the 
Kondo screening cloud have been put forwardi ^d^'^^^d^ 
Also, several studies have emphasized the emergence of 
mesoscopic fluctuations on flnite systems, and the exis- 
tence of even-odd effects in the Kondo cloud when com- 
puted from a lattice model^'^ '^^d^d'^ While there has 
been experimental progress toward the measurement of 
the Kondo cloud )^^i^^ the detection of the spin-spin cor- 
relations has proven to be highly challenging and has 
not been accomplished so far. Depending on the Kondo 
temperature Tk, the Kondo cloud can have a significant 
extension of ~ l^rai^ 

In our work, we examine the spin-spin correlations 
in a real-space model, the single-impurity Anderson 
model (SIAM) that includes charge fluctuations, us- 
ing the density matrix renormalization group method 
(DMRG) i^°i^^i^^ We address two main questions: first, 
we compute the spin-spin correlations between the im- 
purity spin and the conduction electrons at particle-hole 
symmetry and discuss how the Kondo screening length 
can be directly extracted from such data. To that 
end, we discuss several ways of collapsing spin-spin corre- 
lations calculated for different Kondo temperatures onto 



a universal curve. In this analysis, we employ ideas sug- 
gested by Gubernatis et al^ that have also been used in 
previous DMRG studies of the Kondo cloud problemj ^^d^ 
We find that from chains of about L = 500 sites, suit- 
able measures for the L ~ oo screening length can be 
extracted for Kondo temperatures of ksTx/^ ^ 1 • 10^'^ 
(r is the tunneling rate). Knowledge of the universal 
curve further allows us to estimate even for Kondo 
temperatures for which the accessible system sizes are 
too small to host the full Kondo cloud. As a main result 
of our analysis, we find that our measures of S,k extracted 
from the spin-spin correlations have the same functional 
dependence on model parameters as 



(1) 



at particle-hole symmetry (vp is the Fermi velocity in the 
leads, we adopt /cb = 1 throughout the rest of this work). 
The screening length governs the finite-size scaling of 
local quantities such as the polarization or the magnetic 
moment 

Second, wc consider several mechanisms that destroy 
Kondo correlations, namely a gate voltage and a mag- 
netic field applied to the quantum dot. We study the 
changes in the screening length induced by a variation 
in these parameters. We argue that computing the mag- 
netic field dependence of the screening length provides a 
means of extracting the Kondo temperature. 

The emergence of an exponentially small energy scale 
in the Kondo problem, namely Tk, restricts any real- 
space approach with respect to the Kondo temperatures 
that can be accessed. A powerful framework was intro- 
duced by Wilson^^ in the form of the numerical renor- 
malization group (NRG) methodf^iSi which is explicitly 
tailored toward the Kondo problem. This is achieved 
through the introduction of a logarithmic energy dis- 
cretization that allows the Kondo scale to be resolved 
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but loses real space information. Recently, an NRG 
method has been developed to access spatially resolved 
quantities,— li^i^^ extending some older NRG calculations 
for spatially dependent correlation functions Using the 
more recent NRG approach^i^ the spin correlations be- 
tween the impurity and the sites in the leads have been 
computed for the Kondo model, and it has been shown 
that at the Kondo screening length the envelope of 
the correlations crosses over from & decay at dis- 

tances a; < to a l/x decay at distances x > where 
X denotes the distance between the impurity and a site 
in the leads. 

Comparing NRG and DMRG, first, there are technical 
differences between DMRG and NRG with respect to how 
the spin-spin correlations {Si ■ Sj) {Si denotes a spin-1/2 
operator at site i) are obtained. NRG requires a sep- 
arate run for each pair of indices whereas DMRG 
operates directly on real-space leads. That way, after cal- 
culating the ground state for a system of a given length, 
all correlations can be evaluated in a single run. While 
the use of real-space chains is restricted to one dimen- 
sion, which is the case of interest in our work, NRG in 
principle works for higher dimensions too. Second, us- 
ing DMRG, we can gain direct and easy information on 
the finite-size scaling of spin-spin correlations, which we 
heavily exploit in our analysis. Most importantly, DMRG 
can also be applied to quantum-impurity problems with 
interacting leadsi^ that NRG is not designed for. 

DMRG has previously been used to study the Kondo 
cloud in several papers, for both the single-impurity An- 
dcrson modelii and the Kondo model^ii^ In Ref. [ll| by 
Hand et al., in particular, an interesting relation between 
the screening length as extracted from the spin correla- 
tions and the weight of the Kondo resonance has been 
discussed. Our study extends the DMRG literature as we 
consider the mixed-valence regime, the effect of a mag- 
netic field, and we discuss and demonstrate the universal 
scaling of spin-spin correlations for a wide range of pa- 
rameters. Moreover, in the absence of a magnetic field, 
we exploit the SU(2) symmetry of the model in the spin 
sector in the DMRG simulations, which we find is crucial 
for efficiently obtaining reliable numerical results. 

Besides the conceptual interest in understanding the 
scaling properties of the Kondo screening length with 
both system size and Kondo temperature, our results 
are relevant to gauge the range of validity of numerical 
approaches for calculating the conductance of nanostruc- 
tures that employ a real-space representation of the leads 
such as time-dependent DMRG simulations of transport 
in the single-impurity Anderson modeL ^^i^^'^^ Moreover, 
the approaches discussed here to extract the screening 
length could be applicable to more complex geometries 
in a straightforward way, for instance, to multichannel 
and/ or multidot problems. 

Our work is organized as follows. In Sec. |lll we in- 
troduce our model and define the quantities of interest. 
In Sec. mil the spin-spin correlations constituting the 
Kondo cloud arc investigated and we demonstrate how 



to extract the value of the Kondo screening length 
from the spin correlation data, making use of the univer- 
sal finite-size scaling behavior of £,k- Wc proceed with a 
discussion of the behavior of the screening length upon 
driving the system away from the Kondo point via a gate 
potential, presented in Sec. lIVl and then turn to the case 
of a magnetic field in Sec. [V] Wc conclude with a sum- 
mary. Sec. I VII while technical detail on the method and 
computations arc given the Appendix. 

II. MODEL 

We model a quantum dot coupled to a lead by the 
single-impurity Anderson model, describing the lead by 
a tight-binding nonintcracting chain. This constitutes a 
one-channel problem: 

-ff = ^ ednda + BS^ + Uud^Udi 

L-l 
cr 2—1 

Cio- annihilates an electron with spin a =t,i on site i, 
da annihilates an electron with spin a on the dot, and 
n-da ~ dljd^. The spin operators at any site are given 
by Sf = c\^(^1^,c^^, /2, where tr" are the Pauh matrices 
(a = x,y,z). €d denotes the gate potential and B de- 
notes the magnetic field applied to the dot, U denotes 
the strength of the Coulomb interaction on the quantum 
dot, t' denotes the hopping of the dot levels to the first 
site in the lead, t denotes the hopping within the lead. 
The width of the dot level due to the hybridization with 
the lead is given by F ~ 2t'^/t. 

In the absence of a magnetic field, this model has 
a spin SU(2) symmetry. In our analysis, we calculate 
the ground state of this system via DMRG using an 
implcmcntationii exploiting the SU(2) symmetry which 
greatly improves the cfficiency^i^ (see the Appendix 
for more detail). A typical run for L ~ 500 sites with 
m = 1500 states took about 60 h on a 2.6 GHz Opteron 
CPU. 

All simulations, irrespective of e^, are performed at 
half-filling of the full system. As the Kondo scale depends 
exponentially on U /T, while in a real-space representa- 
tion of the leads, the energy resolution is proportional to 
1/L, we restrict our analysis to the intermediate values 
of U /T. The tradeoff for these limitations is that it is 
straightforward to calculate spin correlators, as outlined 
below [see Eq. ^] . 

Throughout this work, we use chains with an overall 
even number of sites. It is well-known that there are 
significant even-odd effects in impurity problem of this 
kindi^iSiiiiii^iii Earlier work (see, e.g.^ Ref. [H), suggests 
that the convergence with system size toward a Kondo 
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Fig. 1: (Color online) Integrated spin-spin correlations T.{x) 
[from Eq. Q] for systems of different sizes at [/ = 1, F = 0.20 
and ed = —U/2. As an example, the threshold of 0.1 that 
we use in Eq. ((SJ to extract ^o.9 is indicated by the dashed 
horizontal line. As an illustration of the typical raw data, we 
show the absolute value of the spin-spin correlations \ {Sd ■ Si)\ 
for L = 300 in the inset. 



state is much faster on chains with an even number of 
sites. We thus work in singlet subspaces. 



III. SPIN-SPIN CORRELATIONS AND KONDO 
SCREENING LENGTH AT ea = -U/2 

In this section, we present our results for the spin- 
spin-correlation function at particle-hole symmetry and 
we discuss two ways of collapsing the data, allowing for a 
determination of the Kondo screening length. In order to 
investigate the behavior of the Kondo screening length, 
we shall study the following integrated spin-correlation 
function, 



{Sd ■ Si) 
{Sd ■ Sd) 



(3) 



to be evaluated in the singlet subspace of the total spin 
Stot = Sd + ^fj^i Si, and under the assumption that 
(5^) 7^ (x is given in units of the lattice constant). 
This definition is motivated by the following convenient 
properties: (i) the decay of S(x) with x characterizes the 
extent to which the total spin of chain sites one to x is 
able to screen the spin on the impurity level, i.e., the 
extent to which X]i=i Si has, crudely speaking, "become 
equal and opposite" to Sd. (ii) When the sum includes 
the entire chain, we always have — 1) = 0. This 
follows by noting that in the subspace with zero total 
spin, where (Sl,) = 0, we have (Sj) = {{Eti S.)^), 
and hence also (S^^) = 2(5^^ • Sd)T.{L - 1). (iii) The 
correlator is normalized to S(0) ~ 1. (iv) In the absence 
of a magnetic field, 'S.{x) is SU(2) invariant, such that 



this symmetry can be exploited in our numerics. In the 
presence of a magnetic field, we shall use a symmetry- 
broken version, replacing {Sd ■ Si)/{Sd ■ Sd) by ((S'JSf ) - 
(5J)(Sf))/((5J5J)-(5J)2). 

As an example, the inset of Fig. [1] shows a DMRG 
result for the absolute value of the bare spin-spin corre- 
lator {Sd ■ Si). The feature at i 200 is a simple effect 
of the open boundary conditions. The spin correlations 
for i smaller than a certain value (here roughly i ~ 200) 
oscillate between negative and positive, while beyond a 
certain point, all {Sd ■ Si) become positive. This feature 
at z ~ 200 precisely appears at the site where this hap- 
pens, i.e., where {Sd ■ Si) with i even changes its sign 
and, as a consequence, the correlator passes arbitrarily 
close through zero. Summing up the correlator according 
to Eq. ([3]) yields S(a;), plotted in the main panel. 

The notion of a screening length is based on the 
premise that the decay of S(x) follows a universal form 
characterized by a single length scale, £^k, as long as this 
scale is significantly shorter than the system size, (,k "C 
L. (According to the expectation that oc Hvp/Tk, 
this condition is equivalent to the following statement: 
perfect spin screening in a system of finite size L can only 
be achieved if the level spacing, which scales like Hvp/L, 
is smaller than Tr-.) Whenever this condition is not met, 
the shape of the decay of S(a;) with x deviates from its 
universal form once x becomes large enough such that 
the finite system size makes itself felt [via the boundary 
condition S(L - 1) = 0]. To extract £,k from DMRG 
data obtained for finite-sized systems, we thus need a 
strategy for dealing with this complication. Below, we 
shall describe two different approaches that accomplish 
this, both involving a scaling analysis. 

To check whether the screening length obtained us- 
ing either of the two scaling strategies conforms to the 
theoretical expectations, we shall check whether its de- 
pendence on the parameters U, T, and agrees with 
that of the length scale = ^ [Eq. Using the 

known form of the Kondo temperature Tk for the An- 
derson modelr^i^ this dependence is given by: 



: exp 



2UT 



(4) 



We shall indeed find a proportionality of the form = 
P^K^ where the numerical prefactor p reflects the fact 
that the definition of Tk involves an arbitrary choice of 
a prefactor on the order of one. Wc emphasize, however, 
that our determination of will be carried out without 
invoking Eq. rather, our results for will turn out 
to confirm Eq. (|4]) a posteriori. In the present section 
we shall focus on the symmetric Anderson model (e^ = 
— U/2) at zero magnetic field, considering more general 
cases inllVI 
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Fig. 2: (Color online) (a) Rescaled integrated spin-spin cor- 
relations E(a;), collapsed onto a universal curve via suitable 
choices of S^k- (b) Comparison of the U and F dependence of 
and [from Eq. ©]: VUT^k (symbols) and VUTpT-^ 
(lines) plotted vs U/T, using p as the fitting parameter (re- 
sulting in p = 6.8). 



A. Scaling collapse of S(a;) 

The first way of extracting the screening length is to 
plot S(x) versus x/^k, where is treated as a fitting 
parameter, to be chosen such that all the curves collapse 
onto the same scaling curve [see Fig. [2]. When attempt- 
ing to collapse the S(a;) data, one faces two issues. First, 
the S(a;) data are nonmonotonic in L, due to the fact 
that the sign oi {Sd ■ Si) oscillates, and for curves scaled 
by different values of £,k, the oscillations are stretched 
by different amounts on a semilog plot. This introduces 
some "noise" to the T,(x) curves, making it somewhat 
difficult to decide when the scaling collapse is optimal. 
Second, for some parameter combinations, the condition 
<C i is not met, and therefore, perfect scaling cannot 
be expected for all the curves. 

These issues can be dealt with by a two-step strat- 
egy: (i) we start with the curves, which collapse the best, 
namely, those with the smallest U/T ratios. These yield 



the smallest S,k values and hence satisfy the condition 
£,K < L required for good scaling well enough such that 
the shape of the universal scaling curve can be established 
unambiguously (to the extent allowed by the aforemen- 
tioned noise), (ii) We then proceed to larger ratios of 
U/T, which yield larger ^/f, and adjust such that a 
good collapse of S(a;) vs. x/^k onto the universal curve 
is achieved in the regime of small x/£,k, where finite-size 
effects are not yet felt. Thus, knowledge of the universal 
scaling curve allows to be extracted even when the 
condition <C i is not fully met. 

The result of such a scaling analysis is shown in 
Fig. [U^a). A universal scaling curve can clearly be dis- 
cerned, with deviations from scaling evident in the curves 
with large U /T, as expected. Moreover, Fig. El^b) shows 
that the results for extracted from S(a;) scaling agree 
rather well with the parameter dependence expected from 
Eq. ^ for p/ {fiVF)-S.K (with a prefactor oip = 6.8), pro- 
vided that U/T > 2. For smaller U/T, no well-defined 
local moment will form and the premise for Eq. ^ no 
longer holds. 



B. Scaling collapse of ^a{L) 

A second strategy for extracting the screening length, 
following Refs. [6lflll. and [13, is to determine the length, 
say , on which the integrated spin-correlation function 
S(x) has dropped by a factor of a of its x = value (for 
instance, a = 0.9 would signify a 90% screening of the 
local spin). Thus, we define 



£,a{L) = min {x; S(a;) < 1 — a} . 



(5) 



The argument of Ca{L) serves as a reminder that this 
length depends on L, since the boundary condition 
I](L— 1) = always enforces perfect screening for x = L. 
However, once the system size becomes sufficiently large 
{L > ^k) to accommodate the full screening cloud, (,a{L) 
approaches a limiting value, to be denoted by [short- 
hand for ^a(oo)] , which may be taken as a measure of the 
true screening length ^k- This is illustrated in the main 
panel of Fig. [T]for a = 0.9: as L increases, the x-valucs, 
where the T,{x) curves cross the threshold 1 — a = 0.1 
(horizontal dashed line), tend to a limiting value. This 
limiting value, reached in Fig.[l]for L > 300, defines ^o.9- 

Figure [3] shows the L dependence of ^0.9 (-^) for several 
values of U /T ranging from 0.4 to 12.5, and system sizes 
up to L = 500. We observe that £,o.g{L,) reaches its limit- 
ing value for small ratios of U /T, which produce ^0.9 val- 
ues smaller than L = 500. For larger values of U/T, how- 
ever, ^0.9 (-^) does not saturate, implying that for these 
parameters, the true screening length is too large to fit 
into the finite system size.— 

Nevertheless, it is possible to extract the true screening 
length in the latter cases as well, by performing a two- 
step finite-size scaling analysis: (i) for those parameters 
U/T for which £,aiL) has already saturated on a finite 
system, we set = ^^(L = 500), and plot ^a{L)/^a vs 
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Fig. 3: (Color online) System size dependence of ^o.g{L) for 
ed = —U/2. Points represent numerical data; lines serve as 
guides to the eyes. 

L/^a- This collapses all such curves onto a universal scal- 
ing curve. For larger U/T, we rescale the £,a{L) curves 
in a similar fashion, but now using as a fit parame- 
ter, chosen such that the rescaled curves collapse onto 
the universal curve determined in step (i). As shown in 
Fig. (Ha) for a = 0.9, this strategy produces an excellent 
scaling collapse for all combinations of U and F studied 
here. 

The above procedure requires the threshold parameter 
a to be fixed arbitrarily. Qualitatively, one needs a large a 
to capture most of the correlations i.e., £,a{L — > oo) ^ £^K^ 
yet a ought not to be too close to one to avoid bound- 
ary effects in the results. Technically, the calculation 
of is much easier the smaller a is, as less correlators 
{Si ■ iSdot) that are of a small numerical value need to 
be computed to high accuracy (see also the discussion in 
the Appendix). For instance, at U/T = 5 and L — 500, 
^0.9 ~ 112 sites, while ^0.75 ~ 29 sites. 

We have carefully analyzed the qualitative dependence 
of our analysis on the threshold a. First, the universal 
scaling behavior in £,a{L)/£,a is seen for a > 0.6. Us- 
ing too small a value for a ignores the long-range be- 
havior of T,{x). Qualitatively, needs to be close to 
the point, where the decay of the envelope of spin-spin 
correlations changes from a power law with 1/x to 1 /x'^ 
(see Fig. 2 in Ref. [11). Second, it turns out that dif- 
ferent choices of a produce values of that differ only 
by a ([/-independent and F-independent) prefactor p{a), 
as illustrated in Fig. IHb) (symbols). In particular, for 
U/T > 2, all follow the same functional dependence 
on the parameters U and F, satisfying the relation 



expected from Eq. ^ (lines in Fig. [5]). It is obvious that 
yields an upper bound to since p{a) > 1 for all 
choices of a. 
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Fig. 4: (Color online) Results of a ^a{L) scaling analysis for 
ed = —U/2. (a) Scaling collapse of Co.9(i)/Co.9 vs I//C0.9, 
obtained by the two-step scaling strategy described in the 
text in Sec. HIT Bl (b) Comparison of the U and F dependence 
of and [from Eq. (|4])] for several values of a: VUT(,a 
(symbols) and VUTp{a) T^^ (lines) plotted vs U/T, using the 
fit parameters p(ffl) shown in the inset (squares). The dotted 
line in the inset indicates the prefactor p = 6.8 obtained from 
the E(a;) scaling analysis of Fig. [^Jb). 



The only exceptions are the data points at U/T = 12.5, 
for which is too large in comparison to L = 500 to yield 
reliable results. The latter are thus excluded when fitting 
the data to determine the best values for p{a), shown 
in the inset of Fig. |4l 

The inset includes the prefactor p = 6.8 (horizontal 
dotted line) obtained in Sec. IIII A| from Fig. [21 via a 
scaling analysis of S(x) (which has the advantage of not 
involving any arbitrarily chosen threshold). Evidently, 
p = 6.8 is rather well matched by p{Q.9) ~ 6.7, imply- 
ing that the two alternative scaling strategies explored 
above, based on T,{x) and £,a{L), yield essentially identi- 
cal screening lengths for a = 0.9. For the remainder of 
this paper, where we consider 7^ —U/2 ot B ^ 0, we 
shall thus determine the screening length by employing 
^o,q{L) scaling, which is somewhat more straightforward 
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to implement than T,{x) scaling. 



IV. GATE POTENTIAL 
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Fig. 5: (Color online) (a) Kondo screening length ^0.9 vs gate 
potential ed/U for several U/T and L — 500. (b) Dot occu- 
pation (ud) vs gate potential. 

We next investigate the behavior of the Kondo screen- 
ing length while sweeping the gate potential applied to 
the dot. Qualitatively, one expects the Kondo tempera- 
ture to increase upon gating the dot away from particle- 
hole symmetry and eventually, as the dot's charge starts 
to deviate substantially from one, the Kondo effect will be 
fully suppressedi^ Consequently, we expect the Kondo 
cloud to shrink upon varying e^. To elucidate this be- 
havior, we focus on values of U/T ^ 5.6 for which 
£,o.g{L ~ 500) yields a good estimate of the true ^k, as 
demonstrated in Sec. IIIII 

Our results for ^0.9 are presented in Fig. [5l^a). In addi- 
tion, and as an illustration, we plot the dot level occupa- 
tion {ud) = {0\ndt + ndi\0) in Fig. [5{b), where |0) is the 
ground state of the system, obtained via DMRG. As we 
shift the dot level away from the particle-hole symmetric 
point at Cd = —U/2 and thus leave the Kondo regime. 



^0.9 falls off rapidly. This is symmetric in the direction 
of the deviation from the Kondo point. In the regime 
-F one would expect Eq. (jj]) to hold roughly. In- 



< 



deed, for = -J7/4, Eq. (g]) still apphes^i^ while for, 
e.g., Cd = this is not the case anymore. The reason 
is that Eq. is only valid in the Kondo regime with 
{ud) ~ 1. From Fig.l^Ib) we see that the dot occupation 
starts to decrease quickly as we increase Cd from — C//2, 
implying that the magnetic moment decreases as well. 
In the mixed-valence regime, Cd ^ — F, ^0.9 measures the 
strength of the spin-spin correlations not originating from 
Kondo physics. 



V. MAGNETIC FIELD 

The application of a magnetic field is known to destroy 
the Kondo effect and its influence on the density of states 
(DOS) and the conductance has been widely studiedi^i^i 
Here, we investigate how the screening cloud collapses as 
the magnetic moment is squeezed by the magnetic field. 
In the presence of a finite magnetic field the total spin 
S is no longer conserved but only is conserved. Thus 
we are left with a U(l) symmetry for 5*^ instead of the 
SU(2) symmetry for S. As a consequence, much more 
computational effort is needed in order to achieve an ac- 
curacy similar to the zero-field case (see the Appendix 
for detail). 

Our results for (i) the screening length ^0.9 = 500) 
and (ii) the magnetic moment of the dot ^ = {{SdY') — 
(S'J)^ are displayed in Figs. [Hlja) andlHlJc), respectively. 
As the magnetic field is increased but still smaller than 
TftT, there are almost no visible effects in ^0.9 (note the 
logarithmic scale in the figure). Once the magnetic field 
B reaches the order of the Kondo temperature Tk, the 
Kondo effect gets suppressed and the extent of the Kondo 
cloud shrinks rapidly. More precisely, a pronounced de- 
cay of the screening length sets in at B ~ 0.5Tx, in agree- 
ment with findings for the field-induced splitting of the 
central peak in the impurity spectral functioui^ Quali- 
tatively, both the screening length and the magnetic mo- 
ment fi exhibit the same behavior. Note that for small 
U /F, charge fluctuations reduce the magnetic moment 
to lie below the value /j, = 1/4 applicable for the Kondo 
model, which presupposes C//F ^ 1. 

To identify the point at which the Kondo effect breaks 
down, we again study the collapse of results from Fig. [S] 
onto a universal curve. This is shown in Fig. [Gfb), and 
as a main result we find: 



eo.9(i?) 
Co.9(S^0) 



(X f{B/TK) 



(7) 



where f{x) describes the universal dependence on B/Tx- 
We note that due to higher numerical effort for calcula- 
tions with a finite magnetic field (as further discussed 
in the Appendix) our numerical results slightly under- 
estimate ^o.9(^) at J7/F > 5, in particular, at small B. 
This, however, has no qualitative influence on the scaling 
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Fig. 6: (Color online) (a) Kondo screening length ^o.g as a 
function of the magnetic field applied to the dot for L — 500. 
In all panels, Tk is given by Tk = fivp/^K with from 
Eq. dH). (b) Scaling collapse of ^o.9(B)/^o.9(-B ^ 0) vs B /Tk 
(c) Magnetic moment fi = ((Sjf ) - (S^f vs. B/Tk. The 
inset shows the rescaled data fj,{B)/fj,{B 0). 



collapse described by Eq. ([7]). We suggest that an anal- 
ysis analogous to the one presented in Fig. [H] could be 
used to extract Tk for models in which the dependence 
of Tk on model parameters is not known. In such an 



analysis. Tk would be the only fitting parameter, since 
^0.9(^7^ ^ 00) can be determined along the lines of 
Sec. mil and one would obtain Tk up to an unknown 
prefactor, which is independent of U/T. 

By rescaling the magnetic moment data to 
ji{B)/fj,{B 0) as shown in the inset of Fig. ^c) 
we again find a universal curve very similar to the 
collapse of ^o.9(^)/Co.9(^ ^ 0) in Fig. E^b). We thus 
confirm that a collapse of local quantities can be used 
to extract Tk, as previously shown using DMRG4 In 
principle, both a scaling analysis of ^0.9 (^) and /i(-B) 
can be used to extract Tk- Using the analysis of the 
screening length data {£,k) offers the possibility of a 
scaling analysis as outlined in Sec. lIIII to reach parameter 
regimes, where a convergence of the data in L has not 
yet been reached. Moreover, the analysis of ^k directly 
unveils the relevant length scales. 



VI. SUMMARY 

In this work, we studied the spin-spin correlations 
in the single-impurity Anderson impurity model using 
a state-of-the-art implementation of the density matrix 
renormalization group method. We first considered the 
particle-hole symmetric point and discussed two ways of 
collapsing the system-size-dependent data onto univer- 
sal scaling curves to extract a measure of the Kondo 
cloud's extension, the screening length £^k, as a func- 
tion of U/r, or Tk, respectively. The first analysis is 
based on a scaling collapse of the integrated correlations, 
while the second one employs a finite-size scaling anal- 
ysis of the distance £,a{L) from the impurity at which a 
certain fraction a of the impurity's magnetic moment is 
screened. .Ja(-^)/Ca(oo) exhibits a universal dependence 
on L/£^a{oo), independently of the parameter U/T. We 
further showed that for an appropriately chosen value 
of the parameter a, both approaches yield quantitatively 
similar estimates of the screening length. Our results for 
£,K, obtained from either of the scaling analyses, nicely 
follow the expected dependence on U /T. 

As DMRG works in real space, the scaling regime could 
only be reached for U/T = 4 and system sizes of L < 
500, but even for larger U/T ^ 6, a collapse onto the 
universal behavior could be achieved. Note that U/T ^ 4 
is the regime, in which time-dependent DMRG is able to 
capture Kondo correlations in real-time simulations of 
transport^l on comparable system sizes, consistent with 
our observations. 

While NRG is better suited to access the regime of 
very small Kondo temperatures Tk, DMRG efficiently 
gives access to the full correlation function {Sd • Si) in a 
single run. As an outlook onto future applications, we 
emphasize that DMRG allows for the calculation of the 
spin-spin correlations in the case of interacting leads^^ or 
out-of-equilibrium, which is challenging if not impossible 
for other numerical approaches with current numerical 
resources. 
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Fig. 7: (Color online) Convergence of ^o.9 vs. the 
residual norm per site Sr/L for ed = —U/2, extracted 
from ground-state calculations using the SU(2) symme- 
try. For each combination of U /V and L, the number of 
states kept increases for data points from right to left as 
m = 200, 400, 600, 800, 1000, 1200, 1500, 2000, 2500, 3000, ex- 
cept for the case fZ/F = 5.6, L — 500, where no point with 
m = 3000 is shown. 



While the first part of our study focused on the 
particle-hole symmetric point where Kondo physics is 
dominant, we have further analyzed how the screening 
cloud is affected (i) by varying the gate voltage and tun- 
ing the system into the mixed-valcncc regime, and (ii) 
by applying a magnetic field at particle-hole symmetry. 
The latter provides an independent measure of the Kondo 
temperature, through the universal dependence of the 
screening length on Tk/B. 

Note added: while finalizing this work, we became 
aware of a related effort on the Kondo cloud, Ref. [s^, us- 
ing the so-called embedded-cluster approximation, slave 
bosons, and NRG. Their analysis is based on calculating 
the local density of states in the leads, as a function of 
the distance from the impurity. 



Acknowledgments 



We gratefully acknowledge fruitful discussions with 
E. Anda, L. Borda, C. Biisser, E. Dagotto, G. B. Mar- 
tins, J. Riera, and E. Vernek. This work was supported 
by DFG (SEE 631, De-730/3-2, SFB-TR12, SPP 1285, 
De-730/4-1). Financial support by the Excellence Clus- 
ter "Nanosystems Initiative Munich (NIM)" is gratefully 
acknowledged. J.v.D and E.H.M. thank the KITP at 
UCSB, where this work was completed, for its hospital- 
ity. This research was supported in part by the National 
Science Foundation under Grant No. NSF PHY05-51164. 



m 

TO 



100 
80 
60 
40 
20 



U/r=3.1, U=1,Ed=-U/2 4( 



B=3-10"'' Tk 
L=500 



■0.75 
^0.8 
Sn.SS 
^0.9 
SO.925 
^ ^0.95 



A 

A 
A 



V ;^c^ 



■)K- -)K 



H H — I h 



10" 



1 0"" 1 0' 

8r/L 




0.19 



0.17 



0.15 



0.13 



Fig. 8: (Color online) Convergence of (,a{B) and the magnetic 
moment jj, vs the residual norm per site Sr/L for a finite 
magnetic field of B /Tk = 3 • 10"^ and ea = -U/2. For 
comparison, we add the B = data, represented by the open 
triangles, from calculations exploiting the SU(2) symmetry 
(m = 1500 states kept). The calculation with a magnetic 
field (symbols connected with lines) uses the U(l) symmetry 
only (m = 100, 200, 400, 600, 800, 1000, 1500 states kept from 
right to left). The results for the dot's magnetic moment fi 
are also included for comparison (solid diamonds). 



APPENDIX A: NUMERICAL DETAIL 

In this appendix we provide details on our numeri- 
cal method. The DMRG calculations presented in this 
work are challenging for two reasons. First, we model 
the conduction band with a chain of length L that pro- 
vides an energy resolution of 1/L, whereas the Kondo 
temperature becomes exponentially small with increas- 
ing U/T [c.f. Eq. (|4|)]. Second, the spin-spin correlators 
are long-ranged quantities making very accurate calcu- 
lations of quantities necessary that are small compared 
to the unit of energy, t. The parameter controlling the 
accuracy of our calculations is the number of states m 
used to approximate the ground state during the DMRG 
sweeps. Typically, we choose m = 1500 (3000 at most) 
for the calculation of the ground state. This results in 
a residual norm per site^^ a measure for the quality of 
the convergence of the calculated ground state towards 
an eigenstate of the Hamiltonian, Sr = {ipo\{H — E)'^\iIjo) , 
on the order ^ = ©(lO"'^). 

Figure [7] illustrates the m dependence of ^0.9 for two 
values of U/T and two values of L at e^j = —U/2, ob- 
tained from simulations using the SU(2) symmetry. The 
larger the ratio U/T and the bigger the system size L, the 
higher the number of states to, needed to be kept to ob- 
tain a well-converged ground state, see Fig. [71 This can 
be understood as follows: higher U/T implies a smaller 
Kondo temperature, i.e., a larger screening length ^0.9 
and longer-ranged spin-spin correlators {S^. ■ Si). A well- 
converged ground state requires these to be evaluated 
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accurately over the entire range i ^ ^0.97 a-nd hence more 
states need to be kept during the DMRG sweeps. For the 
scahng analysis presented in Sec. HII] (see Figs. [2] and [J]), 
we only used data points that are converged with respect 
to the number of states kept. 

In Fig. [51 we illustrate that the convergence with the 
number of states is greatly accelerated whenever the 
SU(2) symmetry can be exploited. We compare this 
preferable case to the calculations with a magnetic field, 
where the SU(2) symmetry is reduced to a U(l) sym- 
metry. In the figure, we use a small magnetic field of 
B/Tk = 3 • 10-3, such that the results for ^^(-B, L = 500) 
coincide with the results for B = 0, previously obtained 
from the SU(2) calculation. For instance, at L = 500 
by keeping m = 1500 states, 5r ~ 3 • 10"'^ is reached 
in the U(l) case as compared to (5r ~ 2 • IQ— * for the 



SU(2) case. For = 1, F = 0.32, we show that this 
residual norm ensures accurate data for up to a = 0.9, 
while for larger a, our U(l) results are well below the 
corresponding SU(2) ones computed with the same m. 

Pragmatically, in the case of broken SU(2) symmetry, 
one may resort to using a smaller threshold a (instead of 
a = 0.9), for which the convergence with m is faster. As 
we have shown in Fig. 21 can be extracted from 
with 0.6 < a < 0.95 up to a nonuniversal prefactor using 
the schemes discussed in Sec. IIIII 

In contrast to the screening length, the calculation of 
the magnetic moment ^, a local quantity, is much better 
behaved. Thus fj, does not suffer much from the slower 
convergence of the U(l) calculation and converges quickly 
to a high precision (displayed as diamonds in Fig. [5]). 
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